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q> ■ abstract 

We estimate the cluster-galaxy cross-correlation function (£ C g), from the APM galaxy 
(— i . and galaxy cluster surveys. We obtain estimates both in real space from the inversion 

' of projected statistics and in redshift space using the galaxy and cluster redshift sam- 

ples. The amplitude of £ cff is found to be almost independent of cluster richness. At 
£ — 1 large separations, r > 5 /i~ 1 Mpc ( h = i/o/100km s _1 , where Hq is the Hubble con- 

\ stant), £, C g has a similar shape to the galaxy-galaxy and cluster-cluster autocorrelation 

functions. £ cg in redshift space can be related to the real space £ cg by convolution with 
an appropriate velocity field model. Here we apply a spherical collapse model, which 
we have tested against N-body simulations, finding that it provides a surprisingly ac- 
curate description of the averaged infall velocity of matter into galaxy clusters. We use 
this model to estimate j3 (f3=^l 6 /b where b is the linear bias parameter) and find that 
it tends to overestimate the true result in simulations by only ~ 10 — 30%. Application 
to the APM results yields (3 = 0.43 with j3 < 0.87 at 95% confidence. This measure is 
OO ! complementary to the estimates made of the density parameter from larger scale bulk 

^\ . flows and from the virialised regions of clusters on smaller scales. We also compare 

the APM £ cg and galaxy autocorrelations directly to the mass correlation and cluster- 
mass correlations in COBE normalised simulations of popular cosmological models 
and derive two independent estimates of the galaxy biasing expected as a function of 
scale. This reveals that both low density and critical density cold dark matter (CDM) 
models require anti-biasing by a factor ~ 2 on scales r < 2 /i _1 Mpc and that the 
Mixed Dark Matter (MDM) model is consistent with a constant biasing factor on all 
scales. The critical density CDM model also suffers from the usual deficit of power on 
large scales (r > 20 /i _1 Mpc). We use the velocity fields predicted from the different 
models to distort the APM real space cross-correlation function. Comparison with the 
APM redshift space £ cg yields an estimate of the value of fi 6 needed in each model. 
We find that only the low O model is fully consistent with observations, with MDM 
marginally excluded at the ~ 2cr level. 

Key words: Galaxies : Clustering ; Large-scale structure of the Universe ; Clusters 
of galaxies; Cosmology. 
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1 INTRODUCTION 

The spatial cluster-galaxy cross-correlation function £ C g(r) 
is defined so that the probability dP of finding a galaxy in 
the volume element dV at a distance r from a cluster is 

dP = n[l + ^ g (r)]dV, (1) 

where n is the mean space density of galaxies. (,c 3 (r) is 
therefore equivalent to the radially averaged overdensity 
profile of galaxies centred on a typical cluster of galaxies. 



The first measurements of galaxy-cluster correlations were 
made by Seldner & Peebles (1977) who measured the an- 
gular cross-correlations of Lick galaxies around Abell clus- 
ters. Their results suggested that £cg(r) was positive and 
significantly different from zero out to large spatial separa- 
tions of r ~ 100 h~ 1 Mpc. A more recent analysis by Lilje 
& Efstathiou (1988) used cluster redshifts to determine the 
cross-correlation between Abell clusters and Lick galaxies 
as a function of metric separation (w(a), see Section 2.2). 
They found no convincing evidence for clustering on scales 
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r > 20 /i _1 Mpc and concluded that some of the signal seen 
by Selder and Peebles was due to artificial gradients in the 
Lick catalogue. However, Lilje and Efstathiou did find some 
evidence for more large-scale power than predicted by the 
'standard' (i.e. Q — 1, h — 0.5, scale-invariant) CDM model. 
On scales smaller than r < 10 fe _1 Mpc, the real-space £ cg {r) 
recovered by inverting the observed form of w(a) is well fit 
by a power law: 



(2) 



with 7 = 2.2 and r = 8.8 h~ 1 Mpc. 

Direct measurements of the spatial cross-correlation 
function from redshift surveys of galaxies and clusters (Dal- 
ton, 1992, Efstathiou 1993, Mo, Peacock and Xia ,1993, 
Moore et al. 1994) have confirmed that £ cs has a similar 
shape to the galaxy-galaxy and cluster-cluster correlation 
functions, but with an amplitude roughly equal to their ge- 
ometric mean. The analyses of Dalton (1992) and Efstathiou 
(1993) were performed using the original APM cluster sam- 
ple of Dalton et al. (1992) and the Stromlo-APM galaxies of 
Loveday et al. (1992). 

In this paper we use compute £ cg using the Stromlo- 
APM galaxy redshift survey and the sample of 364 clusters 
of Dalton et al. (1994b). We investigate its behaviour as a 
function of cluster richness and compare it with the pre- 
dictions of popular comological models. We also calculate 
(,cg(c, 7r), the cross-correlation as a function of separation 
along and perpendicular to the line of sight. The peculiar ve- 
locity field around clusters influences the shape of £ cg (a, n) 
and is expected to depend on f2, and the mass to light ratio 
around clusters, so that we can extract some information 
about the density parameter and galaxy biasing from our 
measurements. We also carry out some specific comparisons 
with N-body simulations (so including non-linear effects) to 
investigate how the biasing of galaxies is expected to vary 
as a function of scale. 



2 CLUSTER-GALAXY CORRELATIONS IN 
THE APM SURVEY 

2.1 The APM data samples 

In this paper we will use three different data samples, the 
APM angular galaxy catalogue, the APM-Stromlo galaxy 
redshift survey, and the APM cluster redshift survey. The 
APM galaxy survey (Maddox et al. 1990a, Maddox et al. 
1990b, Maddox, Efstathiou and Sutherland 1996) consists 
of angular positions and other information, but not red- 
shifts, for over 2 million galaxies with a bj magnitude limit 
of 20.5. The Stromlo-APM redshift survey (Loveday 1990, 
Loveday et al. 1992a, Loveday et al. 1992b) is a survey of 
1787 galaxies randomly sampled at a rate of 1 in 20 from all 
APM galaxies with bj brighter than 17.15. To construct the 
APM cluster sample, an automated procedure was used to 
select clusters from the angular APM survey (Dalton 1992, 
Dalton et al. 1997). Cluster redshifts were then measured, 
and used to construct an original redshift catalogue of 190 
clusters (Dalton et al. 1992) and an extension to 364 clus- 
ters (Dalton et al. 1994). It is the latter cluster catalogue 
which we will use in this paper. In our analyses we will use 
angular clustering and its inversion to obtain real space 3 



dimensional information, using the parent APM galaxy sur- 
vey and the cluster redshift survey (Section 2.2). We will 
also measure the 3 dimensional clustering directly from the 
APM-Stromlo redshift survey and the cluster redshift survey 
(Section 2.3). 



2.2 The projected cross-correlation function 

The simplest statistic which can be used to constrain £ cg {r) 
is the angular cross-correlation function, w cg {6). However, 
the inversion of this quantity to find £ C9 (r) tends to be 
rather unstable. Since we are using the APM cluster redshift 
survey as our cluster sample we can make use of the clus- 
ter redshift information to determine the projected cross- 
correlation function, w cg (a), as defined by Lilje & Efstathiou 
(1988), where a = cz8/Hq is the metric separation of a clus- 
ter with redshift z and a galaxy at angular distance 8 from 
the cluster centre. As our estimator for w cg (a) we use the 
standard estimator, 



w cg (a) 



N ran Ncajcr) 
N gal N C r{<t) 



(3) 



where CG and CR denote cluster-galaxy and cluster- 
random pairs, respectively. We account for the window func- 
tion of the APM Galaxy Survey on a plate-by-plate basis by 
generating a random catalogue for a single plate with 100000 
random points and excising the regions masked from the sur- 
vey as the catalogue is used for each plate in turn. With this 
method, we can use individual galaxy positions from the 
survey data rather than binned cell counts and so we can 
measure w cg (a) accurately at small scales. On larger scales 
we have checked that this method does not introduce large- 
scale power into our determination of w cg (cr) by comparing 
with the results obtained using cell counts for galaxies and a 
single random catalogue for the whole survey (Dalton 1992). 
A similar method has been used more recently in the anal- 
ysis of the Durham-UKST galaxy redshift survey (Ratcliffe 
et al. , 1997). 

The data for w cg (a) are shown in Figure 1 for the clus- 
ter sample cross-correlated with all galaxies to three differ- 
ent magnitude limits. The error bars shown are obtained by 
dividing the APM survey region into four quadrants and de- 
termining the error on the mean from the scatter between 
the four zones. Given the estimate of the depth of the Lick 
catalogue obtained by Maddox et al. (1990a), we would ex- 
pect the points for mn m = 18.5 to correspond to the results 
for w cg (a) obtained by Lilje & Efstathiou (1988). A com- 
parison of Figure 1 to Figure 9c of that paper reveals good 
agreement over the range of a for which w cg can be measured 
reliably. 

The inversion of projected clustering information in the 
form of w cg (a) to the three-dimensional statistic £ C9 (r) in- 
volves a weighted summation of the w cg (a) points (Saun- 
ders, Rowan- Robinson & Lawrence, 1992): 



™Oj+i)-^Oj) 



In 



, +1 + v / CTj 2 +1 - CT y 



(4) 



where r = o». The factor B in Equation 4 accounts for 
the difference in the selection functions of the clusters and 
galaxies and is defined as follows (Lilje & Efstathiou 1988): 
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Figure 1. The projected cluster-galaxy cross-correlation func- 
tion for APM galaxies and APM clusters. Results are shown us- 
ing APM galaxies with different magnitude limits as indicated 
in the figure. The error bars are determined from the scatter in 
the results derived from four nearly equal area zones of the APM 
survey. 
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Here ip is the selection function of the galaxy survey and yi is 
the redshift of cluster i. The selection function tjj was evalu- 
ated using the the luminosity function parameters obtained 
from the Stromlo-APM survey by Loveday et al. (1992). 
Again we show error bars obtained by inverting the w cg (a) 
estimates from four quadrants of the survey for each magni- 
tude range. The data show excellent agreement between the 
three different magnitude limited galaxy samples used, but 
are not well represented by a single power law. 

For the velocity field analysis in Section 4, we will use 
a fit to the real space APM cross-correlation function. We 
choose to fit an arbitrary function which is the sum of an 
exponentially truncated power law and the linear theory cor- 
relation function shape of a scale-invariant CDM model with 
T = Q.h — 0.2 (denoted as £r=o.2(^) below) normalised so 
that a s = iQ 



6*W = (~ 



j3r T) / r 



e vn + a £r=o.2(r). 



(6) 



The parameter combination r c = 11.7, 7 = 2.6, r\ = 0.6, 
(3 = 1.7 and a = 1.8 gives a reasonable fit to the data 
for all magnitude bins and is plotted as a solid line on 
Figure 2. A fit is necessary because the noise level in the 
real space cross-correlation function becomes rather large for 
r > 10 h~ 1 Mpc. The shape of the fit on these large scales 
is motivated by the shape of the redshift space £ cg (Section 
2.2). We also plot on the same figure the power-law fit (given 
by Equation 0) which Lilje and Efstathiou (1988) find is a 
good approximation to the real space cross-correlation func- 
tion derived from Abell clusters and Lick counts. We can 
see that on large scales, there is no evidence for any signifi- 
cant excess of power over this fit. The APM results therefore 
support the conclusions of Lilje and Efstathiou (1988) sum- 
marized in Section 1. Our APM results are in agreement 
with the cross-correlation of APM clusters and Edinburgh- 
Durham Sky Survey galaxies carried out by Merchan et al. 
(1997), although their estimated errors are large. It is useful 
to note that on small scales, r < 2 /i -1 Mpc, our error bars 
on £ C g(r) are very small, so that we will be able to draw 
some interesting conclusions about galaxy biasing on small 
scales from a comparison with theoretical models (Section 
5). 



r (h ^pc) 

Figure 2. The real space cluster-galaxy cross-correlation func- 
tion for APM galaxies and APM clusters. Results are shown for 
different magnitude limits, and the error bars have been calcu- 
lated as described in the caption for Figure 1. For clarity, the 
m Km = 18-5 and m ;im = 20.5 points have been slightly offset in 
the r direction. We also show a fit to the data with form given 
by Equation and the power law (Equation |^) which Lilje & Efs- 
tathiou find is a good fit to the cross-correlation of Abell clusters 
and Lick counts. 



2.3 The redshift space cross-correlation function 

We estimate £c 9 (s), where s represents the separation of 
cluster-galaxy pairs in redshift space, using the APM sam- 
ple of 364 clusters of Dalton et al. (1994b) and the ~ 2000 
galaxies with redshifts from the APM-Stromlo bright galaxy 
redshift survey (Loveday et al. 1992a). The calculation of 
£ cg (s) differs from the evaluation of the cluster-cluster corre- 
lation function (Dalton et al. 1994b) as the galaxy selection 
function falls very steeply with increasing distance. Weights 
must be therefore be applied to the galaxy-cluster pairs to 



* Where erg denotes the rms amplitude of the mass fluctuations 
in spheres of radius 8 h — x Mpc. 
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recover the mimimum variance estimate of £ cg (s) • This op- 
timal weighting (at least on scales for which £ cg (s) < 1) can 
be shown to be (see Efstathiou 1988, Loveday 1990) 



Wij ~ 1/[1 + A%n{ri) Jl 9 {sij)], 



where 



J3 9 ( S ij) 



(7) 



(8) 



Ti is the distance from the observer to galaxy i and Sij the 
separation of galaxy i and cluster j. To use the formula 
we must predict roughly what J^ 9 (s) will be - here we use 
the weighting function resulting from a linear theory CDM 
power spectrum with F = 0.2 and an amplitude twice that of 
£ g9 (s) measured for APM-Stromlo galaxies (Loveday et al. 
1992). After cross-correlating the two catalogues to find all 
galaxy-cluster pairs, we cross-correlate the clusters with a 
catalogue of 100000 random points. This random catalogue 
has the same boundaries and selection function as the galaxy 
sample. We then use the standard estimator to find £ cg (s): 



N ran N CG (s) 

N ga l N CR ( S ) 



(9) 



where Ncg(s) and Ncr(s) are the galaxy-cluster and 
cluster-random pairs in the bin interval centred on s and 
each are weighted using Wij from Equation [?]. 

Results for the 364 TZ — 50 clusters are shown as trian- 
gles in Fig. 3. It can be seen that the there is some curvature 
in the plot, with a definite break at s ~ 30 — 50 /i _1 Mpc. 
The curve crosses £ C9 (s) = 1 at roughly 9 h~ Mpc, which 
is intermediate between the behaviour of the galaxy auto- 
correlation function, £ 9g (s) (so — 5 ft -1 Mpc, in the notation 
of Equation ^, Loveday et al. 1992) and the cluster auto- 
correlation function £ cc (s) (so — 14 /i _1 Mpc, Dalton et al. 
1994). The error bars, calculated from Poisson statisics and 
the number of cluster-galaxy pairs in each bin, are relatively 
small, indicating that £ C9 (s) will be an interesting statistic 
to compare with theoretical models. We can compare the 
results for TZ = 50 clusters with the solid line in this figure, 
which is the fit to the real space cross-correlation function 
(Equation 6). For the moment, we will note that £c 9 (s) is 
marginally higher than £ cg (r) on large scales, boosted by 
streaming motions (Kaiser 1987) and smaller for separations 
less than ~ 4 /t -1 Mpc due to the effect of the cluster veloc- 
ity dispersion. Uncertainties in the cluster redshifts probably 
play a part in depressing the amplitude of £ C9 (s) on small 
scales, as for many cluster we have redshifts for only 2 or 
3 galaxies (each with their own measurement errors). The 
error in the cluster centre of mass velocity could, therefore, 
be as much as a few hundred km s~ . Dalton et al. (1994a) 
have compared results for APM clusters with many (> 10) 
measured galaxy redshifts to the redshifts of the brightest 
galaxy in each cluster. The rms scatter between the two val- 
ues is 512 kms -1 , which should be higher than the error 
on our cluster redshifts as we use > 2 galaxy redshifts per 
cluster. 

Moore et al. (1994) have found that £ cg for IRAS galax- 
ies and Abell clusters is insensitive to cluster richness. We 
have repeated this type of analysis using our APM samples 
and the results are plotted in Figure 3. We have increased 
the lower richness cutoff from TZ = 50 (the full sample) up 
to 1Z = 80, but as the results show there is no detectable 




100 



s (h 1 Mpc) 



Figure 3. The spatial cluster-galaxy cross-correlation function 
for APM-Stromlo galaxies and samples of APM clusters with 
different lower richness limits. All cluster subsamples are drawn 
from the sample of 364 clusters (APM TZ > 50), so that 243, 
113 and 60 clusters are used in the calculation of the curves for 
1Z > 60, 1Z > 70 and 1Z > 80 clusters respectively. The subsam- 
ples have mean intercluster separations of 36 h —1 Mpc (1Z > 60), 
48 h- 1 Mpc ( TZ > 70) and 59 h^Mpc ( TZ > 80). The solid 
line is a fit (with form given by Equation 6) to the real space £ cg 
results for the APM TZ > 50 clusters. 



change in the shape or amplitude of £ C9 (s), given the er- 



rors. Our different samples have space densities n c 



10" 5 /i^Mpc" 3 (TZ = 
9.3 x 10" 6 /i 3 Mpc -3 



(TZ 



3.5 x 
= 60), 



" 6 /i 3 Mpc _3 



50), 2.2 x 10" 5 ft 3 Mpc" 
(TZ = 70) and 4.7 x 10" 
(TZ = 80). We might expect the cross-correlation function of 
richer clusters to have a higher amplitude, at least on small 
scales, as cluster richness should be related to £ cg within the 
cluster selection radius. This was discussed by Seldner and 
Peebles (1977), who found that the effect was smaller than 
expected for Abell clusters. This is probably because the 
distance indicator used by Abell depends on cluster rich- 
ness; for clusters at a given apparent distance, the richer 
objects would actually be further away, thus depressing the 
amplitude of the angular cross-correlation function. These 
problems should not affect the APM cluster sample as it 
has been designed so that cluster richness does not affect 
apparent distance (Dalton et al. 1997). However, the situa- 
tion here is complicated by the fact that we are measuring 
redshift space clustering. Rich clusters will have their clus- 
tering signal smeared out due to their high velocity disper- 
sions, so that the amplitude of clustering will be more heav- 
ily depressed on small scales than for poorer clusters. The 
underlying situation on small scales is therefore not entirely 
obvious. On larger scales (r > 1 /i _1 Mpc) though, our re- 
sults show that the amplitude of £, c g(s) really does depend 
only very weakly on cluster richness. We will show in the 
next section that this is consistent with model predictions. 

We have calculated £ cg for the APM sample in redshift 
space as a function of pair separation along the line of sight 
(n) and perpendicular to the line of sight (a). The effects of 
peculiar velocities, which distort the pair separations in the 
7r direction are evident in the results for our sample plotted 
in Figure 4. We can see elongation present on small scales 
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Figure 4. The APM cluster-galaxy cross-correlation function 
£.cg(&, i") (calculated using the full sample of 364 clusters) shown 
as a function of pair separations perpendicular to the line of 
sight (<r) and along the line of sight (it). Contour levels are at 
£ C9 = 4, 3, 2, 1, 0.8, 0.6, 0.4, 0.2, 0.1, 0., -0.05. The contour level at 
£,cg = 1 is shown by the heavy line; negative contours are plotted 
as dashed lines. £c 9 (o", it) was calculated using 16 bins in a and 
7r in the interval — 40 h~ 1 Mpc. For clarity, before plotting the 
results in this figure they were smoothed using a moving window 
average (3x3 bins). 

which is caused by random galaxy velocities and redshift 
measurement errors. On larger scales, we can see a break 
in the contours around a ~ 6 /i _1 Mpc,7T ~ 10 /i _1 Mpc 
which could be caused by the infall region around the clus- 
ter. This coherent infall should cause compression of the 
contours along the a axis on larger scales (Kaiser 1987, Lilje 
& Efstathiou 1989), but there does not seem to be obvious 
evidence of this in Figure 4. 

The velocity field around clusters should be dependent 
on the mass distribution, and therefore on the value of fi. If 
we have a simple model for how the two main effects present 
in velocity field arise (small scale dispersion and large scale 
infall) we should be able to use distortions in £ C9 (ct, 7r) to 
derive information on and the amplitude of mass fluctu- 
ations. To do this, we need to know £ C9 in real space and 
to have a model which describes the behaviour of galaxy 
velocities around clusters. We apply both of these to our 
£, cg (a,n) data from the APM survey in Section 4.3 below. 
We will first examine the predictions of theoretical models 
using N-body simulations and use them to develop a model 
of the velocity field around clusters. 



3 FROM SIMULATIONS OF 

COSMOLOGICAL MODELS 

3.1 Clusters and galaxies in simulations 

Our predictions of the form of the cluster-galaxy corre- 
lation function in cosmological models come from study- 
ing N-body simulations. We use catalogues of clusters con- 
structed from simulations which have been used in previ- 



ous papers to study the cluster-cluster correlation function 
(Croft & Efstathiou 1994, Dalton et al. 1994b). The sim- 
ulations are of three different spatially flat universes, two 
CDM models and one Mixed Dark Matter (MDM) model. 
One set of simulations is of "standard" CDM (SCDM), with 
fio = l,h = 0.5 and the other is of low density CDM 
(LCDM) with f2o = 0.2, h = 1 and a cosmological constant 
A, where A = 0.8 x 3Hq. The power spectra for the SCDM 
and LCDM models are as given Efstathiou, Bond & White 
(1992). For the MDM model, we used the form given in 
Klypin et al. (1993) with f2 = 1, h = 0.5, and a massive 
neutrino component contributing Q, u = 0.3. We assume scale 
invariant primordial fluctuations for all models. 

Each simulation contains 10 particles in a box of co- 
moving side- length 30000km s _1 and was run using a P 3 M 
A^-body code (Efstathiou et al. 1985). We use 5 realisations 
of each model with different random phases. 

In this Section we use simulated cluster catalogues con- 
structed from the N-body simulations to have the same mean 
separation as the APM sample (d c = 30 7i _1 Mpc). The clus- 
ters are identified from simulations using a friends-of-friends 
algorithm to select candididate centres. The mass enclosed 
within a fixed radius (in this case 0.5 /i -1 Mpc) of the cen- 
tre of mass is computed and the clusters are ordered by 
mass. Finally, a mass limit is applied to generate a cluster 
catalogue of a specified mean space density. The procedure 
is described in more detail in Croft & Efstathiou (1994). 
We note here that as long as the clusters are defined to be 
collapsed objects, the set of objects identified in the simu- 
lations is insensitive to the selection criterion. For example, 
Gaztanaga, Croft & Dalton (1995) identify clusters as high 
peaks in the density field smoothed with small filters and 
recover essentially the same catalogue of clusters (for a va- 
riety of filter sizes) as our percolation algorithm. We can 
also reasonably expect the positions and mass rankings of 
galaxy clusters in our simulations to be insensitive to the 
details of the galaxy formation process, which in the real 
Universe would turn a large agglomeration of dark matter 
into a galaxy cluster. These properties of galaxy clusters 
make them especially useful for elucidating the nature of 
galaxy biasing and galaxy peculiar velocities. In essence they 
consist of a set of fixed reference points around which the 
galaxy overdensity profiles and galaxy velocity profiles can 
be compared directly with observations. 

We do not attempt to carry out any sort of direct iden- 
tification of galaxies in the simulations. Instead, we calcu- 
late the cluster-particle corrrelation function, £ cp . The APM 
£ cg observations can be compared to £ cp to determine what 
sort of galaxy biasing is needed in each model. According 
to linear perturbation theory, the amplitude of as grows in 
proportion to the growth rate D(t) of linear density per- 
turbations (see Peebles 1980, Section 10). The amplitude of 
the two-point correlation function of the mass fluctuations 
thus grows as erf oc D 2 (t) on scales on which linear the- 
ory is applicable. On the other hand, the rich cluster two- 
point correlation function, is almost independent of time (see 
Croft & Efstathiou 1994) because clusters are rare objects 
and are strongly biased compared to the mass distribution. 
We therefore expect the linear theory growth rate of £ C p 
to lie between these two cases, and to be proportional to 
as, as in the high-peak model of Bardeen et al. 1986. This 
means that the choice of output time (and hence normal- 
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Figure 5. The cluster-mass cross-correlation function for simulations of (a) Standard CDM, (b) Low Density CDM and (c) Mixed 
Dark Matter (see text for the parameters of these models). Results arc plotted in real space (solid lines) and in redshift space (dashed 
lines). All models are normalised to be consistent with the amplitude of fluctuations measured by COBE (ignoring tensor modes), so 
that erg = 1.0 for panels (a) and (b) and ag = 0.67 in panel (c). The model curves in panel (c) have been scaled upwards by a biasing 
factor of 1.5 in order to make the model consistent with the amplitude of galaxy fluctuations. Also plotted in each panel (as triangles) 
is the cross-correlation function of APM Stromlo galaxies and the APM cluster sample described in Section 2.3. The error bars on these 
points have been calculated using Poisson statistics. The dotted lines show a fit (Equation 6) to the real space cluster-galaxy correlation 
function, consistent with measurements from the APM survey and also with the cross-correlation of Abell clusters and Lick counts. 



isation of the model) will make some difference to our re- 
sults. The normalisation we have used is consistent with the 
amplitude of fluctuations inferred from COBE microwave 
background temperature anisotropics (Wright et al. 1994) 
so that as = 1.0 for SCDM and LCDM and ag = 0.67 for 
MDM. If we assume that ag for APM galaxies is close to 1.0 
(Gaztanaga 1995) and density evolution on this scale has 
been linear, then £ cp for SCDM and LCDM should be very 
nearly equal to £ C9 . In the case of MDM, when plotting our 
results we merely scale the curve upwards by a constant lin- 
ear biasing factor, 1.0/ ag, . A discussion of more complicated 
biasing, including variations with spatial scale which might 
be required in some models, is deferred to Section 5. 



3.2 Results 

Our results for these three models are shown in Figure 5, 
both in real space (solid lines) and redshift space (dashed 
lines). Concentrating on the real space results first, it can be 
seen that £ C g(r) exhibits a sudden change of slope on small 
scales in all cases. The SCDM and LCDM models have the 
steepest slopes for this part of the curve, which we might ex- 
pect, given that they have the most small scale power. The 
redshift space results show the usual depression on small 
scales and amplification on large scales, both effects being 
largest in the case of the two Q = 1 models, which have the 
largest particle velocities. As far as a comparison with the 
APM £ C g(r) is concerned, we can see that the usual discrep- 
ancy on large scales with SCDM is evident, but that the 
shape on these scales is consistent with LCDM. Indeed for 
r > 3 h~ 1 Mpc LCDM appears to give a good fit in both 



real and redshift space. On smaller scales, MDM appears 
to give a reasonable fit to the shape of both the real space 
results (Equation 6) and those in redshift space. We should 
be cautious about drawing firm conclusions from this, as the 
MDM model has been simulated without including the ther- 
mal velocities of massive neutrinos and this may affect the 
density profiles and internal structure of the clusters. The 
linearly biased MDM model does also seem to have a rather 
high amplitude on intermediate scales, particularly in red- 
shift space (a point that we will return to in Sections 4 and 
5). In any case, it is possible that the efficiency of galaxy 
formation is different near clusters leading to scale depen- 
dent biasing. We will investigate this possibility in Section 
5. 

The richness dependence of £ C j in the models is shown 
in Figure 6, where we plot £ cg (s) (in redshift space) for sim- 
ulated clusters with a similar range of space densities to the 
4 different APM samples shown in Figure 3 . There is a very 
weak richness dependence in the amplitude of of £ C9 , which 
becomes even weaker on large scales (s > 25 /i _1 Mpc). This 
weak dependence of £ cg with cluster richness is compatible 
with the observational results for the APM samples pre- 
sented in Figure 3, which show no significant dependence of 
£ cg with cluster richness. 



4 REDSHIFT SPACE DISTORTIONS AND 
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Figure 6. The cluster-mass cross-correlation function in rcdshift space for clusters with different lower mass limits in simulations of (a) 
Standard CDM, (b) Low density CDM and (c) Mixed Dark Matter. The curves are labelled with the mean intcrclustcr separation of 
each sample. 



4.1 The spherical infall model 

There exists a regime between streaming of galaxies on large- 
scales and the virialised region of clusters which is impor- 
tant in modelling observations of £ cg (er, 7r). As the density 
enhancement within a few ft _1 Mpc of clusters is greater 
than unity, linear theory is not expected to describe the ve- 
locity field accurately. However, if we assume that clusters 
are spherically symmetric, a solution for the non-linear col- 
lapse of the system can be found which is exact before orbit 
crossing takes place. The solution is obtained by treating a 
proto-cluster with a top-hat profile as if it were an isolated 
Friedmann universe with its own value of Qo (see e.g. Regos 
& Geller 1989 for details). Here, we use a good approxima- 
tion to the exact solution due to Yahil (1985) (also used by 
Lilje & Efstathiou 1989), who gives the following expression: 

non—linr \ _ 1 o ' 6 TT ^ ^( r ) I'm"! 

Vinfall (r) - -g^O gor (l + a( r ))0.25 ' ( 10 ) 

where 8(r) is the overdensity inside radius r, 

= ^TT^ -TO = / CepW^ dx. (11) 

Jo 

By way of comparison, the linear theory prediction of the 
infall velocity is 

«iS>.K(r) = -|ng- 6 ff r«(r). (12) 

Galaxy surveys provide with information on the overdensity 
of galaxies inside radius r and not of the mass. Some assump- 
tions about the relationship between galaxies and mass are 
therefore required to infer a value of fto from a measurement 
of Vi n f a ii(r). Here we use the simple linear biasing picture, 
so that J!j P (r) = J^ 3 (r)/b. We therefore parametrise our 
infall model using fi = fig' 6 /&■ As even the non- linear veloc- 
ity model is not expected to hold in the cores of clusters, 
near and in the virialised region, we choose to truncate the 



expression for v™^ 1 ™ (r) with an exponential e -*M/*c«t ; 
with Scut = 50. (See the comparisons with N-body simula- 
tions in the next subsection for a justification of this value 

Of Scut-) 

One further ingredient in our velocity model is the ran- 
dom velocity dispersion about the smooth infalling flow. To 
make things as simple as possible, we assume a velocity dis- 
persion independent of distance from the galaxy to the clus- 
ter, and independent of direction (whether transverse to the 
line between galaxy and cluster, or parallel to it, for exam- 
ple). We also assume that the velocities are drawn from a 
Gaussian distribution. The one-dimensional velocity disper- 
sion will therefore be parametrised by one number, a v . 

4.2 Direct tests of spherical infall on the velocity 
field in simulations. 

In our analysis we need to assume that the velocity field 
predicted from the average cluster density profile is equiv- 
alent to the average cluster velocity field. This is a non- 
trivial assumption and must be tested in some way (in our 
case we will use simulations) before we can make any claims 
for the reliablity of our results. The spherical infall model 
has previously been used in many studies of the infall re- 
gions around individual rich clusters, particularily Virgo 
(e.g. Yahil, Sandage & Tammann 1980, Yahil 1985). How- 
ever, the assumption of spherical symmetry for any individ- 
ual cluster is difficult to justify empirically. In contrast, the 
average cluster profile measured by £ C g(r) possesses spherical 
symmetry by construction. Here we test the validity of our 
dynamical approximations by comparing our predicted ve- 
locity fields directly with the velocity fields measured around 
N-body clusters. We present results both for the mean infall 
velocity as a function of radius and the dispersion about the 
mean infall. 

The average infall velocities of particles, are plotted as 
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Figure 7. The mean infall velocity of particles around simulated clusters as a function of radius. The filled circles in each panel correspond 
to the infall velocities found in N-body simulations of (a) SCDM, (b) LCDM and (c) MDM. The error bars show the error on the mean 
calculated from the scatter in an ensemble of 5 realisations. In each panel we have also plotted the predictions of linear theory (dashed 
lines) and a non-linear spherical infall model (solid lines), as described in the text. The two pairs of dashed and solid lines show the effect 
of including and not including an exponential truncation at high ovcrdensities as described in Section 4.1. 
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Figure 8. The one-dimensional velocity dispersion of particles around simulated clusters as a function of radius. The filled circles in 
each panel correspond to the dispersion in the component of relative velocity along the line from particle to cluster (ic. dispersion about 
the mean streaming motion plotted in Figure 7). The open circles show the one-dimensional dispersion in the transverse direction. We 
plot the same models as in previous diagrams: (a) SCDM, (b) LCDM and (c) MDM. The error bars again show the error on the mean 
calculated from the scatter in an ensemble of 5 realisations. 



a function of distance from the cluster centre in Figure 7 
(filled circles). We also plot the predictions of spherical infall 
(solid lines) and linear theory (dashed lines). The rapid de- 
crease of infall velocity at small r in the simulations is due to 
particles reaching the boundary of the virialised region. We 
have roughly approximated this effect in our velocity mod- 
els by the exponential truncation described in Section 4.1. 
The pairs of dashed and solid curves show the velocity mod- 
els with and without the exponential truncation. Because of 
the somewhat arbitrary nature of this truncation, we will 
restrict our quantitative use of the velocity model to the 



region r > 2.5 fe _1 Mpc. Evidently, the spherical non-linear 
model provides a much better match to the N-body results 
than the linear. Over much of the range of density contrasts 
relevant to our study of the infall region, the non-linear ve- 
locity model underpredicts the infall by roughly constant 
factor, but this is only ~ 30% for LCDM and ~ 10 - 20% 
for SCDM and MDM. At larger radii, r > 20 /i _1 Mpc we 
see that both linear and non-linear models begin to match 
the N-body results even more closely. We therefore expect 
that when we come to estimate j3 from distortions in the 
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cross-correlation function our results will be fairly close to 
the true value. 

The underprediction of the spherical infall model is in 
agreement with results found by Van Haarlem (1992) and 
Diaferio and Geller (1996). As the velocity field in the im- 
mediate vicinity of galaxy clusters has the potential to be 
complicated, there have been many mechanisms suggested 
that might disrupt simple spherical infall. Among these are 
the presence of shear which should speed up collapse (Hoff- 
man 1986, Lilje & Lahav 1991), which does not appear to 
be happening here, and formation of substructure and/or 
ellipticity of the proto-cluster which slows infall (e.g. van 
Haarlem & van de Weygaert 1993). The study of clusters 
in the SCDM scenario carried out by Villumsen & Davis 
(1986) differs from ours in that they analysed the velocity 
field around individual clusters. Here we are dealing with 
the averaged velocity around all clusters. 

Figure 8 is a plot of the 1-dimensional relative velocity 
dispersion as a function of distance from particle to cluster 
in our three cosmological models. The filled symbols show 
the dispersion about the mean infall motion along the line 
from particle to cluster and the open circles show the one- 
dimensional dispersion in the component of relative velocity 
transverse to this line. It can be seen that our assumption 
that cr v is independent of distance and direction turns out to 
be surprisingly good The plots also show that the transverse 
a v is slightly higher than the radial value, in the infall region 
between r ~ 2 h~ Mpc and r ~ 10 h~ Mpc. This might 
be the signature of some sort of "previrialisation" occuring 
(Davis & Peebles 1977, Peebles 1993). 

To illustrate the complexity of the velocity field around 
the simulated clusters, and to see if we can understand why 
some scenarios fit the spherical infall picture slightly better 
than others, we have plotted the particle velocities around a 
sample of individual clusters. In Figures 9 to 11 we show the 
x and y components of the smoothed peculiar velocity field 
(a 2 h~ Mpc Gaussian filter was used) around 9 clusters for 
each model, in the rest frame of the cluster (plotted at the 
centre of each panel) . We have also plotted the density field 
smoothed with the same filter as a grayscale. We can reach 
several conclusions from studying these plots: 

• The magnitude of the velocity field around MDM clusters 
(and to a lesser extent SCDM) is anisotropic - there are sev- 
eral clusters with small arrows on one side and large ones 
on the other. 

• The larger coherence in the velocity field of the MDM 
model compared to the two CDM models is evident. Also 
obvious is the greater magnitude of the velocities in this 
model. 

•There is substucture in the density field for all models and 
the cluster we are interested in at the centre of the panel 
is often part of an elongated system. The MDM model also 
has a noticeably smoother density field on small scales. 

• All models contain some clusters where the main flow pat- 
tern is not centered on the cluster but continues past it, and 
others have more complex and non-radial flows. The pictures 
do not inspire confidence in the idea of simple spherically 
symmetric infall around each cluster. It seems that this pic- 
ture is only likely to apply after some sort of averaging is 
carried out, at least around each cluster (e.g. as in Villum- 
sen & Davis 1986) or around a statistically denned sample of 
clusters, as described in this paper. It is also worth consider- 



ing that these models do have a relatively high amplitude of 
mass fluctuations, being normalized to COBE rather than 
to give the correct abundance of galaxy clusters (White, Ef- 
stathiou & Frenk 1993). It is possible that lower amplitude 
models would conform better to a smooth infall picture. 

Obvious signatures of radial infall, such as overdensity 
"caustics" at the radius of turnaround from the Hubble flow 
(Regos & Geller 1989) have been difficult to find in observa- 
tions of real clusters and in N-body simulations (van Haar- 
lem et al. 1993). This does not seem surprising given the 
complexity evident in these plots. This is consistent with the 
work of Colberg et al. (1997) who have shown that matter 
infalling onto individual clusters in simulations does so from 
specific directions which are correlated with the surround- 
ing large scale structure. It is possible that modelling the 
velocity and density fields as triaxial ellipsoids (van Haar- 
lem & van de Veygaert 1993) may help in the understanding 
of the flow. In any case the accuracy of our spherically aver- 
aged approximation is good enough to permit a quantitative 
analysis of the redshift space distortions. 

4.3 Estimates of f3 and a v 

We are now in the position to make some predictions of how 
the real space £ cg will be distorted as a function of (3 and a v . 
To do this, we shall convolve (,cg(r) with our velocity model, 
so that £ C g in redshift space, £ C g(c, 7r), is given by (Peebles 
1993, Section 20) 

i cg {o,*) = (V^,)- 1 f-°°Jexp(<;- [(v^-'r(r)y)/a v ]f 

x(l + £ C9 M)-exp(-« 2 /^) 2 ] dv, (13) 

where r = y/ (a 2 + n 2 ) and y = yj (r 2 — a 2 ). To model the 
observed estimates of £ C9 (cr, it), we model £ cg (r) in Equa- 
tion 13 with the fitting function of Equation 6 which is con- 
sistent with the real space cross-correlation function of the 
APM sample (Section 2.2). We also test Equation 13 using 
the N-body simulations and empirical estimates of £, cg (r) 
calculated in real space (plotted in Fig. 5). 

To place constraints on f3 and a v we calculate the ex- 
pected £,cg{<y,ir) for a closely spaced grid of values, over the 
range (3 = — 1.5 and a v = — 1500km s _1 . £ C9 (o", yr) i s cal- 
culated at 16 x 16 points, spanning the range of values for 
tv and a plotted in Figure 4. This is done by convolving the 
real space £ cg (r) with our velocity model using Equation 13. 
In Figure 12 we show contour plots for a few sample values 
of (3 and a v . We then find which of our grid of models has the 
smallest \ 2 by comparing with the observed £ C9 (a, 7r). We es- 
timate confidence bounds on /3 and a v from the distribution 
of A\ 2 over the [/?, <j v \ plane. As the infall model does not 
give an accurate prediction for the velocity in the virialised 
region, we have chosen to exclude from the fit the £ C9 (o", 7r) 
point with the smallest value of a and it. Thus we exclude 
any cluster-galaxy pairs with both a and 7v < 2.5 h, _1 Mpc. 

The \ 2 contours are shown in Figure 13, where we have 
tested the SCDM, LCDM and MDM simulations. For the 
test on the N-body models, we have used all (^ 1000) clus- 
ters in each box (and all 10 6 particles), as well as averag- 
ing over 5 realisations. We can see straight away that we 
are predicting values, for (5 at least, that are consistent 
with the results of our direct comparison of infall veloci- 
ties plotted in Figure 7. For the LCDM model, for which 
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Figure 9. The distribution of particles (with their projected velocities) surrounding 9 clusters taken from a simulation of a standard 
CDM universe ( with erg =1.0). The clusters plotted correspond to those ranked by mass as numbers 2,4,8,16,32,64,128,256 and 512 from 
a box of size 300 h~ 1 Mpc. In each panel, the cluster in question is located in the centre, and the velocity field of particles (in the rest 
frame of the cluster) around the cluster has been assigned to a grid using a 2 h~ x Mpc Gaussian kernel. The X and Y components of this 
velocity field in a slice through the centre of the cluster are shown as arrows. The particles have also been assigned to a grid using the 
same kernel and a slice through the resulting density field is also shown on each panel as a grayscale (linear in density). 



the true value is /3 = 0.39 (Lahav et al. 1991), we find 
(3 = 0.521°;°^ (95% confidence limits). For SCDM and MDM 
we measure (3 = 1.201q;?2 an d P = l-Wto.ig, repectively. 
The error bars on (3 given here are computed by marginali- 
sation over all values of a v . The confidence limits on the pair 
of parameters a v and (3 are plotted as contours in Figure 13. 
This diagram shows that given enough clusters we should 
be able to distinguish between high and low (3 models. This 
constraint on j3 from the infall region around clusters (most 
of the signal comes from r < 10 /i -1 Mpc) is complementary 
both to measurements of the central mass of clusters from 
the virial theorem and to larger scale bulk flow observations 
(e.g. Loveday et al. 1996). As the infall method tends to 
systematically overestimate f3 by a small amount, any mea- 
surement which gives a particularly low value, for example 
one inconsistent with high fi models, will be especially in- 
teresting. 

The values of a v we measure using the maxi- 
mum likelihood fits are 370±|5Jkms _1 , 690ig O km s" 1 , and 
570tj?£kms- 1 for LCDM, SCDM and MDM respectively. 
These are again 95% confidence limits obtained after 
marginalising over all values of (3. Comparison with Fig. 8 



reveals that these values are in good agreement with the 
velocity dispersion of particles in the outer parts of clusters. 

We now turn to the determination of (3 and a v from 
the APM sample. The results of the maximum likelihood 
fitting are shown in panel (d) of Figure 13. We can see im- 
mediately that we have much larger statistical errors than in 
our simulations, as there are only 364 clusters and ~ 2000 
galaxies in the APM sample. The value of (3 we obtain is 
0.43, with (3 < 0.87 at 95% confidence. The tail of proba- 
bilities towards rather high values of a v and the best fit a v 
inferred, 660km s _1 (a v < 1070km s _1 at 95% confidence) 
are probably a consequence of uncertainties in the cluster 
redshifts. There is not therefore a good constraint on the 
velocity dispersion in the outer regions of clusters from the 
APM sample. 

5 MODEL-DEPENDENT COMPARISONS: 
BIASING AS A FUNCTION OF SCALE. 

To calculate how biasing of galaxy fluctuations is expected to 
change with scale for each cosmological model, we calculate 
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Figure 10. As for Fig 9 but for an LCDM simulation (with <T8 = 1.0). 



the particle-particle two point correlation function, £ PP (r) as 
well as the cluster-particle cross-correlation function £ cp (r). 
Given the observational results for £ gg (r) and £ cg (r), we can 
define two biasing relations for each cosmological model, 

^W = {»W/UW (14) 

and 

& C9 (r)=£ C9 (r)/£ cp (r). (15) 

If the efficiency of galaxy formation is not changed by prox- 
imity to clusters, then these quantities, b gg (r) and b cg (r), 
will be equal. We use the estimates of £g 9 (s) calculated from 
the APM-Stromlo redshift survey of Loveday et al. (1992) in 
Equation 14. As the observed estimates of £gg and £ C9 arc 
computed in redshift-space, we compare them with redshift- 
space estimates of £ pp and £ cp computed from the N-body 
simulations. The results for the three different models are 
shown in Figure 14. In each panel we show both b gg (s) and 
bcg(s). It is clear that for SCDM, the required bias must 
increase strongly towards larger scales. For LCDM, the de- 
rived b gg is roughly equal to unity and nearly independent of 
scale. There does appear to be a discrepancy between b cg (s) 
and b gg (s) on small scales, but it is likely that errors in 
the cluster redshifts have suppressed the small-scale ampli- 
tude of £, cg (s) leading to an underestimate of b cg . The MDM 
model requires scale-dependent biasing, at least on scales 



< 5 /i _1 Mpc. The dashed lines in the figures show the ex- 
pectations of linear theory and constant biasing for models 
normalised to the COBE measurements according to the lin- 
ear theory value of as - In the case of LCDM, this agrees well 
with the non-linear simulations, although for MDM we can 
see that non-linear evolution appears have made the matter 
fluctuation amplitude higher, and less biasing is necessary 
than linear theory predicts. 

On smaller scales (s < 4 h~ 1 Mpc), we have seen that 
b cg {s) appears to be systematically lower than b gg (s) for 
LCDM and, to a lesser extent, for MDM. As mentioned 
above, errors in cluster redshifts could be causing this ef- 
fect. Furthermore, the relation between b(s) in redshift-space 
and the underlying b(r) depends on the velocity field in each 
model and so mixes physical scales in a complex way. We 
will therefore turn to the real space information that we 
have available, in the form of £cg(r) estimated by inversion of 
the projected APM cluster-galaxy cross-correlation function 
(Section 2.2). Recovery of spatial correlation functions from 
projected and angular statistics can depend sensitively on 
the luminosity function assumed. In Figure 15 we plot b cg (r) 
obtained when three different magnitude limits are applied 
to the angular galaxy catalogue (note the difference in the 
scale of the absicissa compared to Figure 14). On the small 
scales we are considering here, there is not much difference 
between the results and we can see that both the CDM mod- 
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Figure 11. As for Fig 9 but for a MDM simulation (with <r 8 = 0.67). 



els appear to require some anti-biasing. This is in agreement 
with the study made by Jenkins et al. (1998) who derive an 
estimate of bias in models as a function of scale in real space 
from galaxy-galaxy clustering (using Equation [lij). Jenkins 
et al. use a critical density CDM model which has a lower 
amplitude, and so although that model requires b gg (r) > 1 
on all scales, there is a relative antibias on small scales com- 
pared to large. Our MDM models are reasonably consistent 
with a constant biasing factor, consistent with the results 
shown in Figure 5 . Of course as no "hot" particles have 
been used to simulate the neutrinos directly in these mod- 
els, there is some uncertainty in the accuracy of the cluster 
profiles on small scales. 

It should be added that the models which require anti- 
biasing do so on scales which are too small to be probed 
by our non- linear infall model (r < 2 /i _1 Mpc). On larger 
scales, though, redshift distortions of £ cg can give us more 
information about biasing in our models. For example, we 
would expect the bias factors that result from Equations [l4| 
and [l^ to be the same when measured in real and redshift 
space, at least on large scales where linear theory should 
hold. On examination of Figures 14 and 15, this does ap- 
pear to be the case for LCDM, but there is the suggestion 
that the bias factor for both SCDM and MDM is higher in 
real space than redshift space. This can also be seen from 
the plot of £ cg shown directly in Figure 5. This is presumably 



because the coherent peculiar velocities in these two mod- 
els which boost £ cg in redshift space (Kaiser 1987) are too 
high. This is just another way of saying that the amplitudes 
of our models (which are COBE normalised) are probably 
wrong. Lower amplitudes for these two models, such as those 
given by using a normalisation which gives the correct space 
density of galaxy clusters (White, Efstathiou & Frenk 1993, 
Eke et al. 1996) would presumably match the data better. 
Changing the normalization would of course not remedy the 
fact that SCDM does not match the shape of clustering on 
large scales. 

A potentially more powerful way of comparing the mod- 
els and observations involves using the full £ C9 (o", ir) in- 
formation. To do this, we first make the assumption that 
some form of scale-dependent biasing holds so that £ cg (V) in 
the model being tested and the observations are related by 
Equation [l5| We then assume that the velocity field around 
clusters can be predicted from £c P {r) in the model, and that 
the observed form of £ C9 (cr, 7r) has been distorted by this 
velocity field. We use Equation [l3| to model the distortion, 
with the non-linear infall prediction being computed from 
£ C p(V), and used to distort (,c g (r). At this point, we are guar- 
anteed by construction that the amplitude and shape of the 
scale-dependent biasing is correct, but using different values 
of fl and a v in Equation |l3| will result in differing degrees 
of agreement with the observed £ C9 (o", ir). We carry out % 2 
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Figure 12. A smooth fit (Equation 6) to the observed cluster-galaxy cross-correlation function £, cg (r) convolved with a velocity model 
(spherical infall and constant velocity dispersion- see text) to produce a predicted £cg(<7, t) plot. We show results for 3 different values 
of /3 and 3 different values of the velocity dispersion, cr v . Contour levels are the same as in Figure 4 



fitting of (, cg (a,n) exactly as in Section 4.3, except that we 
are constraining fi required in each specific model, rather 
than p for a general scenario with scale-independent biasing. 

The results of this procedure are shown in Figure 16, 
where we can see that all models prefer a relatively low 
value of Q°' e> , which of course is natural for LCDM but not 
for the others. After marginalising over all values of a v , we 
find the following values for ft (with 95% upper limits), 
0.49 (1.28), 0.38 (0.88) and 0.32 (1.06), for SCDM, LCDM 
and MDM respectively. Given the tendency of the spherical 
infall model to slightly overestimate O 0,6 , these results are 
reasonably significant. The a v results for all models are ~ 
700km s _1 . As before, we can attribute this high value to the 
large errors in the redshifts of individual APM clusters which 
should bias the a v results upwards. It is interesting that the 
SCDM model and to a lesser extent the MDM model are 
expected to give a v almost this high without accounting for 
the uncertainties in the velocities. As we know that these 
errors are present (Dalton et al. 1994a), we can regard this 
is additional evidence against these two models. 

In this analysis, we have two parameters, a v and fi ' 6 , 
each of which probes a different part of the velocity field 
around clusters. If the model is correct then the best fit 
values from the observations should match the values of both 
parameters predicted from the model. We have seen that this 
fitting procedure favours the LCDM model over the other 
two we consider here. However, we have been restricted to 
a narrow set of model parameters and it is possible that by 
varying these it might be possible to improve the agreement 
with observations. For example, a lower amplitude SCDM 
model would probably require a higher value of f2 to fit the 
data and so be more consistent. 



This sort of analysis could also be applied to the galaxy- 
galaxy correlation function, although the results would prob- 
ably not be as reliable. The analysis of Section 4 shows that 
non-linear dynamical model of cluster-galaxy peculiar veloc- 
ities works accurately even on small scales. However, it has 
proved difficult to develop a model of galaxy-galaxy veloci- 
ties of comparable accuracy (see e.g. Hatton and Cole 1997). 
Throughout our analysis we have neglected the possibility of 
velocity biasing. Velocity biasing could, if significant, intro- 
duce systematic errors in analyses of redshift-space distor- 
tions. Unfortunately, theoretical predictions for velocity bias 
are so uncertain (e.g. Summers, Davis and Evrard 1995) that 
it is not yet possible to make reliable models of this effect. 



6 DISCUSSION AND CONCLUSIONS 

The analysis presented in this paper differs from other recent 
work on the density and velocity fields around clusters in 
a number of ways. The CNOC cluster redshift survey of 
Carlberg et al. (1997) specifically targets galaxies within and 
close to the virialised region of clusters. Carlberg et al. find 
that a rescaling of cluster profiles allows all clusters to be fit 
by the universal dimensionless halo profile of Navarro, Frenk 
and White (1997). We do not carry out such a rescaling 
here, as we are interested in the profiles of clusters out to 
a much larger distance (where presumably rescaling would 
not work, as all profiles must converge to the mean density) . 
Our dynamical measure of the dark matter in each cluster 
is also based on a different velocity and density regime, the 
infall region, rather than the virialised region. 

The sort of analysis we employ here could be extended 
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Figure 13. Contours of constant A\ 2 corresponding to 68%, 95% and 99.7% confidence limits on values of a (the 1 dimensional velocity 
dispersion) and /3 (~ Q°' 6 /b) found by fitting a velocity model to £cg(&, it). In panels (a) to (c) wc show the results found from clusters 
identified in 5 ensembles each of simulations of the SCDM, LCDM and MDM models respectively. In panel (d) we show results for the 
APM survey (using the £ C g(cr, 7r) plotted in Figure 4). In each panel the values of a and j3 for which the \ 2 was evaluated are shown by 
dots, and the best fit values are indicated by a cross. 
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Figure 14. The bias b(s) as a function of scale in rcdshift space for 3 cosmological models. b(s) is found by dividing the APM ^cg(s) by 
^cp(s) for each model, the results being represented by open circles in each panel. A separate estimate of b(s) can be found by evaluating 



£gg( s )/£pp an d ttu s quantity is plotted as filled circles. The dotted line shows in each case the relation that would be expected if there 
exists a linear biasing relation between galaxies and mass and <rs for galaxies = 1.0. 
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Figure 15. The bias b cg (r) as a function of scale in real space for 3 cosmological models for the region within 10 ft _1 Mpc of cluster 
centres. b(r) is found by dividing the APM £cg(r) ) by ^cp(r). The real space cross-correlation function is found by inverting the projected 
cross-correlation function of APM galaxies and clusters (see Section 2.2). Results for three different galaxy magnitude limits are shown, 
m Hm = 18-5 (filled circles), mu m = 19.5 (open circles) and m; im = 20.5 (open squares). The amplitude of the models has been adjusted 
to our usual COBE normalisation so that erg = 1.0 for panels (a) and (b) and erg = 0.67 for panel (c). The dotted line shows in each 
case the result for b cg (r) that would be expected if there exists a linear biasing relation between galaxies and mass and if erg for galaxies 
= 1.0. 




Figure 16. A% 2 distributions for Q 6 and a v evaluated in a similar way to those in Figure 13, but this time using the velocity field 
predicted from the mass distribution in models to distort the real space £ cg . This will give an estimate of the values of f2 - 6 and <r v 
needed in the models to account for the redshift distortions, given that they are automatically correctly biased to give the same real 
space £cg as the APM survey. 



to constrain O from the virialised region. In the usual anal- 
ysis of cluster velocity dispersions, there is no way of telling 
unambiguously whether a given galaxy is a cluster member 
with a high velocity or a foreground or background galaxy. 
Non-member galaxies and the various ways employed to 
prune them can add significant bias to determinations of 
the velocity dispersions from samples of clusters (see e.g. 
van Haarlem et al. 1997). Analysis of the "finger of God" 
effect on the cross-correlation function offers a way of treat- 



ing the problem statistically, which automatically accounts 
for contamination. In this paper, we have not made use of 
the information on <r v because we have a relatively small 
number of galaxies in our sample and because there are un- 
certainties in the errors of the APM cluster redshifts. Future 
surveys such as the 2dF and Sloan redshift surveys (Colless 
1997, Gunn & Weinberg 1995) and also the presently avail- 
able Las Campanas Redshift survey (Schectman et al. 1996) 
could be used effectively for this purpose. These large sur- 
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veys will also form good datasets to which we can apply the 
non-linear infall models used in this paper, and should pro- 
vide accurate estimates of Q and biasing parameters from 
the infall region of clusters. 

The value of /3 we obtain in this paper is consistent with 
that obtained from a consideration of linear redshift distor- 
tions of the galaxy-galaxy clustering function on larger scales 
(e.g. Loveday et al. 1996) . This suggests that any significant 
biasing as a function of scale is probably confined to scales 
smaller than ~ 2 /i _1 Mpc. 

In summary, we have investigated the shape and ampli- 
tude of cluster profiles both in the APM survey and in sim- 
ulations of cosmological models with Gaussian initial condi- 
tions. In all cases we find a distinct two-component struc- 
ture, with the shape of £ cg (r) on large scales consistent with 
the galaxy-galaxy or particle-particle correlation functions. 
For r < 5 /i _1 Mpc we see a steepening of the profile. This 
region is also the part most sensitive to differences in the 
richnesses of clusters used to calculate t; cg {r), although the 
effects are too small to be noticeable in the observational 
sample. 

As expected for objects forming in a bottom-up sce- 
nario, the velocity and density field around clusters in the 
simulations is complex, with substructure evident in plots. 
Nevertheless, when we average the infall velocities around 
all clusters, the shape and amplitude of the infall curve is 
surprisingly well described by a spherically symmetric col- 
lapse model. It is therefore feasible to use this simple veloc- 
ity model to describe distortions in (, cg (<J, 7r) and therefore to 
constrain Q. When we do this using the APM data, we find 
that f2°' 6 /b has a best fit value of 0.43 and is constrained 
to be less that 0.87 at 95% confidence, where b is the linear 
bias parameter. 

We have also compared £ 93 and £ cg with N-body sim- 
ulations of the mass distribution to estimate the biasing 
of galaxy fluctuations as a function of scale. We conclude 
that anti-biasing on scales r < 2 /i -1 Mpc is required in the 
SCDM and LCDM senarios, as their cluster profiles appear 
to be steeper than in the observations. This sort of biasing 
might be possible, as galaxies or quasars forming in clus- 
ters might suppress galaxy formation nearby (e.g. Babul & 
White 1991). On the other hand, this might be hard to rec- 
oncile with a galaxy density field which is totally unbiased 
on large scales (as may he required in LCDM). The MDM 



model goes have the correct cluster profiles, without any 



scale-dependent biasing. However, we find that the redshift 
distortions in the cross-corrrelation function predicted by 
this model are too strong, by about 2a, so that the LCDM 
model should still be considered the most successful of those 
tested. 

Application of the techniques presented in this paper to 
larger redshift surveys should yield robust measures of the 
cosmic density from the infall region around galaxy clus- 
ters. The spherical infall model is a dynamical approxima- 
tion which is accurate in the highly non-linear regime and 
should continue being useful for constraining Q and biasing 
in a fashion complementary to traditional linear analyses of 
galaxy-galaxy clustering. 
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